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(g) Procede et dispositif permettant a un modem de se synchroniser sur un transmetteur de donnees 
numeriques par voie hertzienne en presence de brouilleurs. 



(57) Le proc6de consiste a estimer (9) la matrice 
de correlation ftxx des signaux recus du trans- 
metteur sur un ensemble de n capteurs de 
reception, a calculer (10) la matrice de correla- 
tion inverse a calculer (11) des vecteurs 
d'intercorrelation f^j entre les signaux recus 
sur I'ensemble des n capteurs et un signal de 
replique d(k) connu, a calculer un critere de 
synchronisation (15) en effectuant les produits 
des vecteurs d'intercorrelation transposes 
conjugues ftii et de la matrice des vecteurs 
d'intercorrelation fxdi. et & comparer (16) la 
valeur du critere obtenu a une valeur de seuil r\ 
detenminee pour placer la synchronisation sur 
rechantillon pour lequel la valeur des criteres 
depasse la valeur de seuil in* Pour estimer le 
decalage Doppler, la sequence de synchronisa- 
tion est 1 decoupee en tranches d'un nombre 
determine de symboles pour effectuer des cal- 
culs de transformation de Fourier rapide sur 
des vecteurs d'intercorrelation determines sur 
chaque tranche. 
Application : Transmission radio de donnees. 
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La presente invention cencerne un proc6de et un dispositif permettant a un modem de se synchroniser 
sur un transmetteur de donnees numeriques par voie hertzienne notamment en presence de brouilleurs. 

Dans les systemes de transmissions de donn6es numeriques par voie hertzienne la synchronisation est 
effective, k partir de sequences de symboles connus du recepteur et intercaiees entre les symboles d'infor- 
5 nation, par correlation entre le signal recu et un signal de reference construit k partir des symboles connus. Le 
resultat est compare k un seuil pour decider de la synchronisation. Des systemes permettant la synchronisation 
d'un modem k partir d'un seul capteur et d'un seul recepteur sont par exemple connus de rarticle D.V. Sarwate et 
M.B. Pursely intitule "Cross correlation properties of pseudo random and related sequences" paru dans IEEE Vol. 
66, pp. 593-619. May 1980. Mais les techniques mises en oeuvre qui permettent en outre une estimation des ca- 
10 racteristiques du canal de transmission deviennent rapidement inefficaces en presence de brouillage. 

Lorsque le rapport signal k bruit n'est pas suffisant pour la prise de synchronisation par une technique 
mono-capteur, des techniques visant k accroTtre le rapport signal £ bruit et a permettre la synchronisation sont 
utilisees. 

Ces techniques peuvent etre regroup6es en deux families. Une premiere famiile vise a pr6traiter le signal 
15 re?u avant d'effectuer la synchronisation du signal pretraite. La seconde famiile vise k ameliorer le rapport 
signal k bruit et k effectuer la synchronisation conjointement. 

Dans les techniques relevant de la premiere famiile il faut egalement distinguer deux categories : 
: ^ V . une premiere categorie met en oeuvre un pretraitement mono-capteur de type excision de brouilleurs 
' ? avaht d'effectuer la synchronisation. Ce traitement est connu par exemple de rarticle de L.B. MILSTEIN 

20 "Interference rejection Techniques in spread spectrum communications", paru dans Proc: IEEE, Vol. 76. 

pp. 657-661. June 1988. Suivant cet article on cherche k pr6dire la contribution du brouilleur. dont la 
bande est suppos6e plus petite que la baride du signal, sans pr£dire celle du signal utile. Cela permet 
d'ameiiorer les performances de la synchronisation en presence d'interf6rences bande etroite (type 
FSK). Mais le traitement reste ineff icace en presence d'un brouillage occupant toute la bande utilisee 
25 par le signal utile. Cela est en partie due au fait que les techniques de pretraitement mono-capteur ne 

permettent un traitement des signaux que dans un espace limite k I'axe des temps et k I'axe des fre- 
quences. 

- La seconde categorie concerne les techniques mettant en oeuvre un pretraitement multi-capteurs de 
rejection, precedent la prise de synchronisation. L'idee consiste a mettre en ceuvre un filtrage spatial. 
30 discriminant le signal et les brouilleurs par leur direction d'arrivee, le temps, I e spectre ou la puissance 

comme cela est decrit dans Touvrage de R.A. Monzingo, T.W. Miller, intitule "Introduction to adaptive 
arrays", edite par John Wiley and Sons, New-York. 1980 et la these de P. Chevalier, intituiee "Antenne 
adaptative : d'une structure lineaire a une structure non Iin6aire de Volterra". These de doctorat, Uni- 
versite de Paris Sud, Juin 1991 puis de correier le signal de reference avec la sortie de I'antenne. Les 
35 techniques de filtrage spatial discriminant le signal et les brouilleurs par le temps, le spectre ou la puis- 

sance sont generalement sous-optimales, alors que celles exploitant la connaissance de la direction 
d'arriv6e du signal necessitent la mise en oeuvre d'une technique de goniom6trie, couteuse et delicate 
dans la gamme des frequences HF notamment du fait des multi-trajets. ■ 
Pour ces raisons, les techniques de la seconde famiile, visant a effectuer la rejection de brouilleurs conjoin- 
40 tement avec la prise de synchronisation, a partir de la connaissance du signal de reference, apparaissent plus 
prometteuses. Une telle technique a par exemple ete proposee par Compton dans son article "An adaptive 
array in a spread spectrum communication system" paru dans Proc. IEEE, vol. 66, pp. 289-298, March 1978, 
pour la synchronisation de signaux PN, k partir d'une architecture analogique integrant des "Delay Lock Loop". 
Une seconde technique fut propos6e par L.E. Brennan et I.S. Reed dans leur article intitule "An adaptive 
45 array signal processing algorithm for communications", paru dans IEEE Trans-Aero-Syst. vol. AES-18, n°1, 
pp. 1 24-1 30, June 1 982 pour les transmissions numeriques et compietee par D.M. Duglos et R.A. Scholtz dans 
rarticle intitule "Acquisition of spread spectrum signals by an adaptive array", paru dans IEEE Trans-Acous- 
Speech-Signal-Proc, vol. ASSP-37, n°8, pp. 1253-1270, August 1989. Ces techniques apparaissent toutefois 
beaucoup plus couteuses que les techniques mono-capteiir et sont malheureusement inapplicables telles quel- 
50 les en presence des derives Doppler du signal inhere ntes aux transmissions ionospheriques. 
Le but de invention est de pallier les inconv6nients precites. 

A cet effet, I'invention a pour objet un proc6d6 permettant k un modem de se synchroniser sur un trans- 
metteur de signaux numeriques par voie hertzienne en presence de brouilleurs consistant k: 

- estimer la matrice de correlation Rxx de signaux regus sur un ensemble de N capteurs de reception par 
55 blocs composes d'un nombre d'echantillons d'une sequence de synchronisation transmise par le trans- 
metteur, , . A 

- k calculer la matrice de correlation inverse R^J, 
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- a calculer des vecteurs d'intercorreiation f Xd entre les signaux X(k) recus sur Tensemble des N capteurs 
et un signal de replique d(k) connu, 

- a calculer un critere de synchronisation effectuant le produit scalaire entre le yecteur d'intercorrelation 
r X( j et le vecteur obtenu par produit de I'inverse de la matrice de correlation Rxx 1 et du vecteur d'inter- 
correlation r Xd 

- a comparer la valeur du critere obtenu a une valeur de seuil r\ determinee pour placer la synchronisation 
sur I'echantillon pour lequel la valeur des criteres d6passe la valeur de seuil ti, caracteris6 en ce qu'il 
consiste pour estimerle decalage Dop pier in herentaux communications entre emetteurs etrecepteurs, 
a decouper la sequence de synchronisation en un nombre determine de M tranches de ko symboles, a 

effectuer sur ces tranches des caiculs de vecteurs d'intercorrelation partielles r Xd i 0 = 1 M) de telle 

sorte que sur chacune des tranches de ko symboles, le decalage Doppler fasse tourner la phase d'une 
valeur inferieure a 90°, a effectuer des caiculs de transform6e de Fourier rapide (FFT) sur les vecteurs 
fxdi obtenus, et a calculer une valeur du critere de synchronisation a partir des vecteurs r Xd (kAf), obtenus 
en resultat des caiculs de transformee de Fourier rapide, et de la matrice de correlation R^. 

L'invention a egalement pour objet un dispositif pour la mise en oeuvre du procede. 

Le procede et le dispositif selon l'invention ont pour avantages qu'ils permettent a la fois d'assurer la syn- 
chronisation d'un modem dans des conditions de brquillage severe et en presence de multi-trajets, d'estimer 
la derive eh frequence entre emeUeur et 'r^cepteiir et de detecter les instants d'arrivee de tpus les multi-trajets 
associ6s au signal utile, pourvu que ceux-ci soient de puissance suff is amment forte par rapport a la puissance 
de bruit de fond. 

D'autres caracteristiques et avantages de l'invention apparattront dans la description qui suit faite au re- 
gard des figures des dessins annexes qui represented : 

- La figure 1 un schema fonctionnel d'une antenne adaptative selon l'invention. 

- La figure 2 un exemple de mise en oeuvre d'un f iltre adapte spatial a discrimination par la modulation 
pour la realisation d'une antenne adaptative selon l'invention. 

- La figure 3 un schema d'antenne adaptative pour la mise en oeuvre du procede de synchronisation spa- 
tiale selon l'invention. 

- La figure 4 un graphe pour illustrer la methode d'estimation de la matrice de correlation de signaux pro- 
venant de N capteurs. 

- La f igure 5 un graphe pour illustrer le decoupage des symboles de la sequence de synchronisation par 
tranche pour le calcul de vecteur d'intercorrelation en presence de decalage Doppler. 

- La figure 6 un graphe montrant revolution de critere de synchronisation mis en oeuvre par ('invention 
en presence d'un decalage Doppler. , 

- La figure 7 un organigramme retracant I'ensemble des etapes necessaires a la mise en oeuvre du pro- 
cede de synchronisation selon l'invention. , 

Le procede de synchronisation mis en oeuvre par ('invention fait appel aux techniques connues de syn- 
chronisation mono-capteur pour les appliquer k la realisation d'antennes adaptatives sur replique. 

Selon la technique de synchronisation mono-capteur, le critere de synchronisation C(n) est obtenu en cal- 
culant la correlation y(n) entreK symboles connu s dans une sequence de synchronisation et les K echantillons 
qui suivent I'arrivee d'un echantillon n surle capteur, en normalisant les r6sultats par la puissance p(n) du signal 
recu. Ces caiculs sont effectues en appliquant les relations : 

y( n > = ^S x (k + n)d*(k) (1) 
K k=i . . 

i - r - ■ < , 

P( n ) = -^ £x(k+-n)x*(k + n) (2) 
k=i 

.... ci "'='S |! ,3) 

- x(k) est le signal issu du capteur utilise, 

- d(k) est le signal de r6 pi ique connu (K symbole$ places en debut de trame), 



# 
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- * represente I'operation de conjugaison sur les nombres complexes. 

Ces calculs sont reiteres pour chaque echantillon n, et lorsque le critere C(n) depasse un certain seuil ?i, 
dependant de la probability de fausse alarme, la synchronisation est decidee. L'instant de synchronisation est 
determine parle rang d'echantilldn n 0 pour lequel C(no)> n- La normalisation par la puissance p(n) permet de 
choisir un seuil independant du niveau de signal re$u. 

Dans ce precede la synchronisation est obtenue lorsque le rapport signal £ bruit est suff isant. Ses per- 
formances dependent en partie du nombre d'echantillons pris en compte pour la correlation. Cependant en 
presence de brouilleurs forts, la synchronisation ne peut etre effectuee. 

L'idee mise en oeuvre dans Hnvention consiste alors a effectuer la synchronisation, en sortie d'un traite- 
ment d'antenne adaptative, permettant la rejection des brouilleurs. 

Le schema fonctionnel d'une antenne adaptative est represente k la figure 1 ,- celle-ci comporte un reseau 
1 de N capteurs 1 ! ... 1 N . un reseau de f iltres adaptes 2 (constitute par des circuits de ponderations amplitude- 
phase, des f litres £ reponse impulsionnelle f inie RIF ou inf inie Rll) couple Sun processeur de commahde adap- 
tatif 3; Le processeur 3 calcule les coefficients Wi a W N du reseau de f iltres 2 pour I'adapter au signal utile 
et aux brouilleurs, en I'absence de toute connaissance a priori sur les brouilleurs et h partir d'une information 
minimale sur le signal utile. Les coefficients \N, a W N permettent de former un trou dans le diagramme de 
rayonnernent du reseau d'antenne dans la direction des brouilleurs, tout en maintenant un certain gain dans 
la direction du signal utile. " , * : ; ■ ;; 

4 - Differerits types d' information sur le signal utile sont expioites par I'antenne adaptative. Ceux-ci sontd'or- 
dre spatial lorsque la direction d'arrivee du signal utile est connue, d'brdre temporelle lorsque les instants de 
presence et d'absence de signal utile sont connus par exemple en evasion de frequence, reiatifs a la forme 
d'onde lorsqu'une information sur la modulation du signal utile est disponible pour disposer d'une replique du 
signal utile en intercalant par exempie des sequences connues dans la trame utile, d'ordre spectral lorsque 
le signal utile occupe une certaine bande de frequence et que les brouilleurs occupent une bande plus impor- 
tante, ou statistiques pour separer les differentes sources arrivant sur le reseau en exploitant par exemple les 
statistiques d'ordre 4 des signaux en supposant que les differents signaux sont statistiquement independents. 

Dans le cas de I'antenne adaptative representee & la figure 2, 1'information sur le signal utile est constitute 
par les sequences de symboles connus intercalees dans la trame entre les sequences de symboles d'infor- 
mation. Les f iltres places en sortie de chacun des capteurs 1 , & 1 N sont formes respectivement par des circuits 
de ponderation amplitude-phase W, & W N des signaux X1 (t) a x n (t) fournis par les capteurs. Les signaux obtenus 
en sortie des circuits de ponderation sont appliques a un circuit de sommation 4. 

Les caracteristiques du reseau de f iltres 2 sont calcuiees en minimisant une Erreur Quadratique Moyenne 
(EQM) estimee £ entre un signal de replique d(k) et le signal y(k) sortant de I'antenne de la fapon schematisee 
& la figure 2 ou les elements homologues & ceux de la figure 1 sont reperes avec les memes references, § est 
def inie par la relation 



K K ' ' ' 

i - Z|y(K) - atop = z w l x(k) - d(k)f (4) 

k=1 k=1 



ou: 

- y(k) est la sortie de I'antenne adaptative et est def inie par 

y(k) = WX(k) (5) 

- W est le vecteur poids def inissant I'antenne adaptative, W 6tant le vecteur transpose conjugue corres- 
pondant. (le signe t en exposant def init Toperateur correspondant.) 

- X(k) est le vecteur forme par les signaux repus sur les capteurs. 

- d(k) est le signal de replique connu. 

La difference y(k) - d(k) est obtenue sur la figure 2 & la sortie d'un circuit soustracteur 5. 
Le vecteur poids W est calcule par la formule connue de Wiener : 

W = R-jAd -(6)- 



avec : 



RxX=^i;X(k).X*(k). (7) 
K k=1 



45 



55 
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matrice de correlation des signaux provenant des capteurs 1, a 1 N . (X l (k) etant le vecteur transpose conjugue 
du vecteurX(k)). 



? Xd = ^£x(k)d» 
K k=i 



vecteur d'intercorrelation entre X(t) et d(t) (8). 

10 L'antenne adaptative ainsi utilisee conduit a une maximisation du rapport signal/(bruit + interference) 

connu sous la designation SNIR abreviation anglo-saxonne de "Signal to Noise + Interference Ratio" obtenu 
en sortie de l'antenne, avec une possibility de rejection pouvant aller jusqu'a N-1 brouilleurs, et est appelee 
Filtre Adapte Spatial (FAS). 

La recherche de la synchronisation est effectuee selon ('invention en sortie de l'antenne adaptative en 

15 discriminant le signal selon sa modulation, et en utilisant le critere de synchronisation mono-capteur. Le sche- 
ma de principe de la synchronisation spatiale est represents a la figure 3 ou les elements homologues a ceux 
des figures 1 et 2 sont representes avec les m ernes references. Le critere de synchronisation, devient, en uti- 
lisant les notations precedentes : 

20 

C(n) = l|y(k + n)d(k) (9) 
K k=i 

25 y(k+n)represente la sortie de l'antenne adaptative conduisant a une minimisation de I'erreur quadratique 

moyenne 4(n) estimee sur K echantillons entre la replique d(k) et le signal arrivant pris a partir de I'echantillon 
n, suivant la relation 

K ' ' 

30 £(n) = £|y(k + n) - d(k)p avec y(k + n) = W(n)\ X(k + n) (1 0) 

k=1 

avec': ' ■ • A 

.35 W(n) = Rxx- 1 (n)-f X d(n) (11) 

Les differentes correlations sont calculees sur les K echantillons qui suivent I'arrivee d'un echantillon n 
suivant les relations : 

Rxx(n) = ^ix(k + n)X l (k + n) (12) 



1 K 

rxd(n)=-prX x ( k + n) d *( k ) 03) 
K k=i 



Lorsque I'echantillon n correspond a la position de synchronisation, les poids W(n) def inissent le filtre 
so adapte, et y(k+n) contient alors le signal utile debarrasse de la contribution des brouilleurs. Dans le cas contrai- 
re, aucun signal n'est correle avec la replique d(k) et y(k+n) ne contient plus alors que la contribution du bruit 
de fond. 

L'expression du critere de synchronisation spatiale se simplif ie en ecrivant : 



K ( K 

C(n) = -1 XWlnl'xfk^nJdfk^Wln) 1 1 £X(k + n) d*(k) (14) 
K k=1 ' l K k=1 



5 
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soit A A A 

C(n) = ^ Xd (n)Rxx(n)- 1 rxd(n) (15) 

De meme que pour le cas mono-capteur, ce calcul est reitere pour cheque 6chantillon n, et lorsque te critere 
C(n) depasse un certain seuil dependant de la probability de fausse alarme. la synchronisation est d6cidte 
(la normalisation par la puissance est ici obtenue par I' intermediate de R«). 

L'instant de synchronisation est alors determine par le rang de l^chantillon pourtequel C(n 0 ) est superieur 
a une valeur de seuil determinee r\. 

Les calculs correspondent sont effectues sur la figure 3 par le bloc 7 de calcul du critere C(n) et par le 
d^tecteur de seuil 8, £ partir de la sequence de synchronisation d(t) et du signal y(t) corr6les par le correlates 
6 

Cela n6cessite pour chaque 6chantillon. d'une part de calcuier N correlations entre le signal de replique 
d(k) et les signaux x(k) recus par les N capteurs (soit N fois les calculs effectu6s pour la synchronisation mono- 
capteur). et d'autre part, d'estimer puis d'inverserla matrice d'intercorrelation des signaux capteurs. 

Au cours de ces calculs la synchronisation et le calcul du f iltre adapts 2 sont effectues conjointement. 
Lorsque la valeur du critere dSpasse le seuil r\, calcute, la synchronisation est decid6e et le vecteur poids d£- 
finissant le f iltre 2 s'en dSduit facilement par la formule connue de Wiener. Le vecteur poids W peut ensuite 
etre utilise pour filtrer les 6chantillons contenant Tinformation. ■ 

. En supposant que la replique d(t) est de puissance 1, les cqrr&ations tendent vers les valeurs suivantes, 
lorsque lenombre d'ecriantillons K devieht trfes graiid : 1 • - 

fxd — VttT S 
Rxx - R + n s S S* 

ou: 

- R est la matrice de correlation des interferences et du bruit de fond, 

- S est la matrice du vecteur directeur du signal utile, 

- k s est la puissance du signal utile recu surun capteur. t 

Le critere de synchronisation qui converge alors vers une valeur C 0 s'Scrit encore/ 

r = *» S>R:1 S = SNIRopt (16) 
Co 1 + ^ StR-i S 1 + SNIRopt 
oil SNIRopt est le SNIR asymptotique obtenu en sortie du f iltre spatial adapts 2. 

Lorsque le SNIR est fort, le critere tend vers 1, lorsqu'il estfaible, tl tend vers 0. L'importance du choix du 
seuil n. fixe la probability de detection de la synchronisation. Le seuil depend naturellement de la probability 
de fausse alarme souhaitee. A 

Le calcul du critere de synchronisation C f de f Xd et de Rxx est effect pour chaque position de la syn- 
chronisation et notamment la matrice Rxx est estimSe f et invers6e egaiement pour chaque position de la syn- 
chronisation, la puissance de calcul demandee est 6lev6e. Mais il est possible de r6duire la puissance de calcul 
en estimant la matrice de correlation Rxx moins souvent sur un intervalle de temps plus important (K' 6chan- 
tillons au lieu de K avec K'>K de la facon representee a la figure 4) et en utilisant cette estimation de Rxx sur 
toutes les positions de synchronisation incluses dans cet intervalle de temps. 

Dans ces conditions la matrice Rxx s'6crit : 



Rxx K . 



k=1 k=K 0 



+ e. (17) 



ou 



- B(t k ) est la contribution du bruit et des brouilleurs sur X(k), 

- S{t k ) est la contribution du signal utile sur Xfo), 

- Kq est Tinstant d'apparition du signal utile sur le palier, 

- e est un terme prenant en compte I'intercorrelation entre le signal B(t) et le signal utile S(t), qui tend vers 
0 lorsque K devient grand: 

Pour la position de synchronisation testee : ' 



6 
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fXd = K 



K 0 4K 
k=Ko 



(18) 



10 



15 



20 



25 



Dans les conditions asymptotiques (avec K»1), lorsque le signal utile est synchronise avec la replique 
d(t), Rxx et f Xd se simplif ient et s'ecrivent : 

Rxx^R + K '-^- + 1 n s SSt (19) 

fxd^V^S (20) 

Le critere obtenu en estimant la matrice de correlation Rxx sur K' echantillons possede comme valeur 
asymptotique, apres un calcul utilisant le lemme d'inversion matricielle: 



C 0 = 



SNIRopt 



+ *"-*> + 1 SNIRopt 



(21) 



K' 



K' — Ko +1 
Comme - — < 1 , alors C' 0 > C 0 



La valeur du critere apparait dans ces conditions augmem d'apps- 
rition du signal utile sur la frequence consideree, ce qui entrame une amelioration de la prise de synchronisa- 
tion, au prix d'une puissance de calcul reduite. 

D'autre part, le nombre K' ne doit pas etre choisi trop grand, sa valeur est determinee pour que sur 
K' echantillons, I'environnement de brouillage puisse §tre considere comme stationnaire. 

Le seuil a prendre en consideration depend de la probability de fausse alarme (Pfa). La (Pfa) est la pro- 
bability qu'un bruit soit detecte comme une synchronisation. Le calcul de la valeur du seuil n a lieu en calculant 
Pfa=P(C>r)) lorsque I'antenne recoit uniquement du bruit. 

En presence de bruit, le critere de synchronisation s'ecrit: 



30 



35 



c = 




(22) 



En supposant que B k est un bruit gaussien complexe de puissance a 2 , independant entre les capteurs, 
lorsque K est suffisamment grand, le terme du milieu de la relation (22) devient : 



40 



^ K k=1 



= a +2 Id (23) 



45 



(Id etant la matrice identite ou unitaire) 

le critere de synchronisation peut alors etre approxime de la fa9on suivante : 



K 2 a 2 



( K 

u=1 



\ l ( K 



(24) 



50 



Soit : 



55 



N 



C = 



w 2 2 



k=1 



2^ 



K CT i=1 



avec: 
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ZBj^^Rj-hjl; (26) 
k=1 

Bt represente le bruit obtenu sur ie capteur i, pour I'echantillon k. 

R t + jlj peut etre approxim6 par un bruit gaussien complexe de variance K a 2 . Dans ces conditions Rj et l { 
10 sont chacun des bruits gaussiens r6eis de variance 



15 



En notant T, = Ri ou jj, C devient: : .. 



20 



a 2N, v 

02 ^3 2(f) < 27 > 



25 



K o 2 

ou T, est une variable aleatoire qui suit une loi normale centre de variance — -— . La variable aleatoire Z f = 



V2 



T| suit done une loi normale centree de variance 1 . 



C devient : 



30 



- 2N 

C^SZ? (28) 



i=1 



35 



Comme la suite des Zj peut etre consid6r6e comme une suite de 2N variables al6atoires independantes 
gaussiennes de moyennes nulles et de variances 1, la variable aleatoire 



40 



2N 
i=1 



45 



suit une loi du % 2 a 2N degres de liberty. 
La probability d'avoir Csn s'6crit alors ; 



P(C > Tl) = P 



50 



f2N - 

V>1- 



= Pfa (29) 



Le tableau suivant donne differentes valeurs de seuils r) obtenues par recherche dans une table donnant 
P(ZSA) pour une loi du x 2 a 2N degres de liberie et pour Pfa = 1 0" 5 . 
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15 



N 
N 
N 
N 



5 
4 
3 
2 



A = 42 
A = 37 
A = 33 
A = 28 



ti = 0.26 
ti = 0,23 
ti = 0,21 
11 = 0,18 



Comme entre remission et la reception, le signal transmis peut etre affecte d'un certain decalage en fre- 
quence, ou decalage Doppler, du essentiellement aux derives en frequences des oscillateurs locaux des 
synthetiseurs a remission et a la reception, ainsi qu'aux conditions de propagation ionospherique, il est ne- 
cessaire de compenser ce decalage. Le risque est s'il n'est pas compense, d'empecher toute prise de syn- 
chronisation. 

En appelant Af 0 , le decalage Doppler ehtre le signal recu par I'antenne et le signal attendu et en considerant 
la position de synchronisation, le vecteur signal recu peut s'ecrire avec les notations des paragraphes sui- 
vants : 

X(t) = ^ d(t) e 2 i*-*V S + B(t) 
La matrice de correlation estimee peut alprs s'6crire r:. - • :■ J~y- »(..■■■.■ . 



20 



RXX = ^ Z [V^s" d(t k )e 2 i" Af o t K S + B(tk)J^s" d(t k )e 2 -i' t Af o «K.S + B(t k )] (30) 
k —1 
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35 



40 



45 



RxX = e + 



ou e est une matrice prenant en compte les intercorrelations entre la contribution du signal utile et le bruit b(t). 
30 Lorsque K»1, dans les conditions asyrnptotiques, Rxx peut s'approximer par : 

Rxx = n s SS + + cr 2 ld (32) 
La relation (32) montre que le decalage Doppler ne modifie pas le calcul de En effet, le decalage 
s'elimine naturellement lorsque le signal utile est multiplie par sa quantite conjuguee. 
Le vecteur d'intercorrelation peut s'ecrire : 



? Xd = £ Z[^ d(tK)e 2 ^ S + B(t k )] cf (t k ) (33) 



? xd - V + 1 i[V^Kt k )^ a^*-** s] (34) 



oil e" est un terme prenant en compte I'intercorrelation entre le bruit B(t) et la replique. Lorsque K»1 , r X(J peut 
s'approximer par : 



50 



k=1 



(35) 



La relation (34) montre que le decalage Doppler modifie I'estimation du vecteur d'intercorrelation r Xd et 
que du fait que la phase du terme en exponentiel complexe peut tourner de plusieurs fois 360° sur I'intervalle 
de temps utilise pour le calcul de la correlation, Integration ne peut plus avoir lieu de facon coherente ce qui 
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empeche le critere de synchronisation multicapteurs de converger. Le processus suivant permet d'estimer le 
decalageDoppleretdelecompenser pour rendre possible la prise desynchronsation. 

Selon ce processus, les K symboles de la sequence de synchronisation sontdecoupes de la facon repre- 
sentee a la figure 5 en M tranches de Ko symboles (M K 0 = K). de telle sorte que sur chacune des tranches 
de Ko symboles, le decalage Doppler fasse tourner la phase du signal recu par I'antenne d une valeur infeneure 
a 90° af in de permettre une integration coherente de r X(J sur chacune des tranches. 

Si r x<ll designe la valeur de i'intercorrelation estimee sur la tranche i, celle-ci s'ecnt : 
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r Xdi = 



(1+1)K 0 

^ k=iK 0 +1 



Af 0 



^|d(kT)f e F * 



S i = 0,. ...M-1 (36) 
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En mettant en facteur dans la relation (36) un terme dependant uniquement de Ko, la longueur de la tranche 
elemental consideree, et dir decalage Doppler, ia reiat^ deyient : 



20 



r Xdi = 



2j7t.^.iK 0 

b F « .S i = 0,...,M-1 (37) 



25 



30 



Ce qui peut encore s'ecrire 



fxdi = x, e^- 1K o .S i = 0,.,.,M - 1 (38) 



La relation (38) montre que les vecteurs r Xdi sont des vecteurs, echantillonnes a ^. Leur fransformee de 
Fourier discrete est alors : 



35 



40 



45 



M-1 



Af 



WAf) s Z r m -e ' Fe/Ko ' 
i=0 



i=0 



M-1 2jz Af °- Af 



Fe 



iKo 



(39) 



(40) 



La relation (40) montre que lorsque Af = Af 0 , le decalage Doppler est compense. ce qui rend maximum la 
valeur du critere de synchronisation r Xd (Af) t Rxx" 1 rx<j(Af)- /. * , , • a ™ 

Af in d'estimer Af. c'est-a-dire af in d'estimer le decalage Doppler qui conduit a la.valeur maximale du cn- 
tere et afin de reduire le cout de calcul correspondant, une transformee de Fourier rapide est effectuee sur 
les vecteurs rxd. (le nombre de tranches elementaires M etant choisi pour former une puissance de 2) : ce calcul 
donne : ' 



50 



M-1 



-2jn, 



Af, 



-.ik M-1 



ivi— i — .tjii.— — — .1* ■»•-■ 

fcdO^sZ r Xdi e = £ *di e 

i=0 i=0 



-2 m - 

"M k = 0 M-1 (41) 



oil Afi est la resolution Doppler telle que : 

Af 1= ^ .(42) 

Si la resolution est insuf f isante, une amelioration d'un facteur 2 peut etre obtenue en calculant une trans- 
formee de Fourier rapide, sur les vecteurs r xd , "decalees" en frequence de la facon su,vante : 



10 



EP 0 665 665 A1 



r Xd kAf, + 



\ M-1 
J i=0 



-j7l.— 

r Xdi e M 



_ 2jn ik 
e *M k = 0 M-1 



(43) 
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^operation consiste a multiplier chacun des vecteurs d'intercorrelation r Xd i par le terme e- d"™) et a effec- 
tuer ensuite un calcul de transformee rapide sur les vecteurs ainsi obtenus. 

Le critere de synchronisation est ensuite estime sur chacune des cases de calcul des transformees de 
Fourier precedentes par la relation : 



15 



20 



c[kf)=fi d (kf)R^f xd (kf) k = 0 2M-1 (44) 

Bien evidemment, si la resolution est toujoursinsuffisante, il est possible de I'ameliorer en utilisant le 
meme processus pour calcul er des vecteurs /-'A-:.//.. '■■ ' 
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ou n est un entier quelconque. 

Un test sur la case de la transformee correspondant au maximum du critere de synchronisation est effec- 

tue. Si le seuil est depasse pour une valeur k<> ^ donnee du decalage Doppler, il est alors possible d'aff iner 

le calcul du decalage Doppler (indispensable pour la demodulation) en calculant la valeur du critere pour dif- 

ferentes valeurs de decalage Doppler entourant ko En effet, a Tissue de cette etape, la resolution sur I'es- 

timation du decalage Doppler est Pour I'ameliorer, il faut proceder par dichotomie : a chaque iteration et 

de la facon representee a la figure 6, la valeur du critere pour la valeur de decalage Doppler situee au milieu 
des deux valeurs de decalage Doppler donnant les valeurs maxima du critere lors de ('iteration precedente 
est calculee. La recherche par dichotomie est arretee des que la resolution souhaitee sur I'estimation du de- 
calage Doppler est obtenue. Apres n iterations de Talgorithme de recherche par dichotomie, la resolution, qui 

etait initialement de est : . 



Af t / 2 = Afj 
2 n 2 n + 1 



(45) 



L'enchamement des etapes de calcul precedente pour la mise en oeuvre du procede selon I'invention a 
I'aide de moyens de traitement appropries formes par exemple par des processeurs de traitement du signal 
ou tout dispositif equivalent connu est figure sur I'organigramme de la figure 7. 

II consiste suivant Tetape 9 a estimer la matrice de correlation Rxx suivant ia relation (7). Cette matrice 
est estimee sur un bloc contenant un nombre d'echantillons K ? superieur au nombre d'echantillons de la se- 
quence de synchronisation. 

L'estimation de Rxx est utilisee ensuite pour toutes les positions de synchronisation incluses dans le bloc. 
Puis la matrice de correlation inverse R^ 1 est cateulee a I'etape 10. f Xdi qui estdefini par la relation (8) est cal- 
cule & I'etape 11. Pour chaque echantillon correspondant a une position de synchronisation P dans la suite 
des echantillons P = 1 a P = K'-K, f Xd i est calcule pour i = 1 a M sur M blocs de Ko symboles comme defini a 
la figure 5 suivant la relation 
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fxdi = ^- S ■ X(P + K + iK 0 ) d (K + iK 0 ) (46) 
K o K=1 
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Des calculs de transformee de Fourier Rapide sont ef tectums £ I'etape suivante 12 sur uri nombre deter- 
mine M (M = 8 par exemple) de vecteurs f Xdj suivant la relation (41 ). 

Aux etapes 13 et 14 un calcul de vecteur R Xd i decale et un calcul de transformee de Fourier rapide suivant 
la relation (43) sont executes pour ameliorer la resolution par un facteur de 2. 

Le critere de synchronisation est calcule a I'etape 1 5 en application de la relation (44). 

La comparaison du entire ou seuil ti est effectuee a I'etape 16. Si k I'etape 16 la valeur du critere 




est superieure ou egale au seuil, la position de synchronisation est definie a I'etape 17 sur la position P de 
I'echantillon correspondent L'estimation du decaiage Doppler peut ensuite etre aff ine & I'etape 1 8 en calculant 
h nouveau la valeur du critere de synchronisation pour des decalages Doppler entourant le decaiage 

^^ f ...,.-=- : ^ l : : ^,. ,., ..... . 

precedent. Cette etape peut etre suivie par une etape 19 d'estimatioh de la position de synchronisation pour 
les autres trajets. 

La recherche de synchronisation precedemment decrite est effectuee pour chacun des echantillons de 
signal numerise. En pratique, pour ameliorer la precision de la position de la synchronisation il est preferable 
d'effectuer un surechantillonnage du signal d'un facteur determine L, L etant egal & 4 par exemple. Les trai- 
tements de I'organigramme de la figure 7 et decrits ci-dessus sont encore utilisables sur le signal ainsi sure- 
chantillonne, mais il est possible d'utiliser une variante dans le procede, af in de reduire la puissance de calcul 
necessaire au traitement. Dans ces conditions, apres avoir calcule pour chaque echantillon le vecteur r Xdi pour 
i = 1, ... M sur M blocs de Ko symboles suivant la relation (46) un critere de synchronisation C, est calcule sur 
chacun des M blocs par la relation 

C| = ^-R^rxdi (47) 
Puis un critere somme de I'ensemble des C| est calcule suivant la relation : 



Ci(p) = X c i ( 48 ) 

i=i; 

Ce critere correspond a une integration non coherente des Ci (cbntrairement & ce qui est fait pour le calcul 
par FFT ou on tntegre de fagon coherente), ce qui se traduit par une remontee du critere sur les echantillons 

de bruit. , 

Le calcul du critere a lieu par blocs de L echantillons consecutifs. Sur chacun de ces blocs, I'echantiUon 
qui donne la valeur maximum de C^P) est selectionne pour effectuer le calcul par transformee de Fourier ra- 
pide necessaire & l'estimation du decaiage Doppler comme dej& decrit aux etapes 113 19 de la figure 7. 



Revendications 

1. Procede permettant £ un modem de se synchroniser sur un transmetteurde signaux numeriques par voie 
hertzienne en presence de brouilieurs consistant & : 

- estimer (9) la matrice de correlation Rxx de signaux recus sur un ensemble de N capteurs de recep- 
tion (1i ... 1 N ) par blocs composes d'un nombre d'echantillpns d'une sequence de synchronisation 
transmise par le transmetteur,. A 

- a catculer (10) la matrice de correlation inverse R^, 

- a calculer (11) des vecteurs d'intercorrelation f Xd entre les signaux X(k) recus sur I'ensemble des N 
capteurs et un signal de replique d(k) connu, ^. % . . 

- a calculer un critere de synchronisation^) ef fectuant le produit scalaire entre le vecteur d'intercorre- 
lation r Xd et le vecteur obtenu par produit de I'inverse de la matrice de correlation R^ et du vecteur 
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d'intercorrelation r X( j 

- a comparer (8) la valeur du critere obtenu a une valeur de seuil ti d6termin6e pour placer la synchro- 
nisation sur l'6chantillon pour lequel ia valeur du critere depasse la valeur de seuil t|, caracterise en 
ce qu'il consiste pour estimer le decalage Doppler inherent aux communications entre emetteurs et 
recepteurs, a decouper la sequence de synchronisation en un nombre determine de M tranches de 
k 0 symboles, a effectuer (11,13) sur ces tranches des calculs de vecteurs d'intercorrelation partielles 
f Xd i (i = 1 . M) d e telle sorte °. ue sur chacune des tranches de ko symboles, le decalage Doppler 
fasse tourner la phase d'une valeur inferieure a 90°, a effectuer (12, 14) des calculs de transformee 
de Fourier rapide (FFT) sur les vecteurs f Xd i obtenus, et a calculer une valeur du critere de synchro- 
nisation a partir des vecteurs r^kAf^ obtenus en resultat des calculs de transformee de Fourier ra- 
pide, et de la matrice de correlation Rxx- 

A 

Precede selon la revendication 1 caracterise en ce que ^estimation de la matrice de correlation Rxx est 
obtenu par la relation 

■■ ■■■■ •• ' k=l 
et le vecteur d'intercorrelation est estime suivant ia relation 



r Xd = 7rZX (k+p) d[ k)> 

k=1 

ou p est la position du premier echantillon a partir duquel la correlation est effectuee entre le vecteur X 
et la replique d. 

Precede selon la revendication 2 caracterise en ce qu'il consiste a choisir la case de calcul de la trans- 
formee rapide qui donne une valeur maximum du critere de synchronisation et a comparer la valeur maxi- 
mum obtenue a une valeur de seuil determinee pour decider de la synchronisation. 

Precede selon Tune quelconque des revindications 1 a 3 caracterise en ce qu'il consiste ( pour obtenir une 
meilleure resolution sur le decalage Doppler et une valeur plus importante du maximum du critere, a cal- 
culer (1 3) des vecteurs r xdi decaies en frequence et a effectuer un calcul de FFT sur les vecteurs r xdi de- 
cays de facon a obtenir la valeur du critere avec une meilleure resolution sur les frequences Doppler. 

Proc6de selon Tune quelconque des revendications 1 a 4 caracterise en ce qu'il consiste a surechantil- 
lonner les signaux recus sur les N capteurs par rapport au rythme symbole, pour ameiiorer la resolution 
sur la position de synchronisation. 

Precede selon I'une quelconque des revendications 1 a 5 caracteris6e en ce qu'il consiste a effectuer le 
calcul de la matrice de correlation Rxx sur un nombre d'echantillons K' sup6rieur au nombre K d'echan- 
tillons de la sequence de synchronisation af in de diminuer la puissance de calcul. 

Precede selon I'une quelconque des revendications 1 a 6 caracterise en ce qu'il consiste pour diminuer 
la puissance de calcul et apres avoir surechantillonne d'un facteur L par rapport au rythme symbole, a 
effectuerune etape de pretraitement consistant a calculer pour chaque echantillon p (p = 1 a L) composant 
un symbole les valeurs C ( du critere de synchronisation sur chacune des M tranches de ko symboles et 
a calculer la valeur du critere en effectuant la somme des criteres C, obtenus pour retenir dans la suite 
du traitement I'echantillon donnant la valeur maximale du critere C = ZC,. 

Dispbsitif pour la mise en oeuvre du precede selon I'une quelconque des revendications 1 a 7 du type 
comportant, un r6seau de N capteurs (1 v „1 N ) de signaux de reception couple a un reseau de f iltres adap- 
t6s (2) comportant des circuits :, de ponderations amplitude phase commandes par un processeur de 
" cdmmande (3), le processeur de comm'ande (3) etant programme pour 
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estimer la matrice de correlation inverse des signaux re$u X (K) et calculer les vecteurs d'inter- 
correlation f xd entre les signaux X(k> repus etle signal de replique d (K) caracterise en ce qu'il comprend 
des moyens de traitement (5.. .19) pour calculer un critere de synchronisation (7) en effectuant un 
produit scalaire entre le vecteur d*intercorr6lation r xd et le vecteur obtenu par le produit de I'inverse 
de matrice de correlation et du vecteur d'intercorrelation rxd. 

comparer (8) la valeur du entire obtenu & la valeur de seuil t| pour placer la synchronisation sur 
l'6chantillon pour lequel la valeur du critere d6passe la valeur du seuil t|, 

et estimer le d6calage Doppler inherent aux communications entre emetteurs et recepteurs, en d£- 
coupant la sequence de synchronisation en un nombre determine de M tranches de ko symboles, 
en effectuant (11, 13) sur ces tranches des calculs de vecteurs d'intercorrelation partielles f xdi (i = 

1 M) de telle sorte que sur chacune des tranches de Ko symboles le decalage Doppler fasse tour- 

ner a la phase d'une valeur inf6rieure d 90°, en effectuant (12, 14) des caiculs de transform6e de 
Fourier rapide (FFT) sur les vecteurs f X d, obtenus, et en calculant une valeur du critdre de synchro- 
nisation & partir des vecteurs r X d(KA^, obtenus en r6sultat des calcuis de transformee de Fourier ra- 
pide et de la matrice de correlation Rxx- 
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Calcul du critere : 
C(n)=-J-z: y(k+n)d"{k] 

C ln)=r Xd {n)R" xx (n)r Xd (n) 



Synchronisation si : 

C(n) >n 

(r\ dependant de la Probability 
de fausse alarme desiree) 
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Processeur 
W(n) = Ri(n)r Xd (n) 
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FIG.3 
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Calcul du critere C (k^ 0 ) 
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